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Abstract 

We consider a two-particle system in which a particle is held fixed, and the 
other one rotates around the axis perpendicular to the line joining the parti- 
cles centers. The rotating particle leads to a displacement of its polarization 
charge on the surface. Our results show that the rotational motion of the 
particles generally reduces the force between the particles. The dependence 
of interparticle force on the angular velocity of rotation will be discussed. 
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I. INTRODUCTION 



The prediction of the strength of the electrorheological (ER) effect is the main concern 
in the theoretical investigation of ER fluids [1—5]. An ER fluid is a suspension of particles 
with high dielectric constant (conductivity) in a host medium of low dielectric constant 
(conductivity). The ER effect originates from the induced interaction between the polarized 
particles in an ER fluid. Upon the application of an external electric field, polarization 
charges are induced on the particle surfaces, leading to anisotropic forces between the par- 
ticles. Under the influence of the ER effect, the particles in an ER fluid initially aggregate 
rapidly into chains within milliseconds and then the chains aggregate into columns within 
seconds [3,4]. The rapid field-induced transition between the fluid and solid phase makes this 
material both important for wide industrial applications such as shock absorbers, dampers 
and clutches, as well as for experimental and theoretical investigations. 

In deriving the induced forces between the particles, existing theories assume that the 
particles are at rest [6-10]. In a realistic situation, the fluid flow exerts force and torque 
on the particles, setting the particles in both translational and rotational motions. For 
instance, the shear flow in an ER suspension exerts a toque on the particles, which leads 
to the rotation of particles about their centers [11]. Recent experiments showed that the 
induced forces between the rotating particles are markedly different from the values predicted 
by the existing theories which have not taken the motion of particles into account [12]. 

Wang designed a delicate experiment, in which a spherical particle is held fixed and the 
other one rotates uniformly, to investigate the dynamic ER effects ||13| . The experimental 
result is that the force exerted on the rotating particle is smaller than that on the same 
particle at rest. In order to explain this effect, we formulate a model, which describes the 
relaxation process of the polarized charge on the surfaces of the rotating particles. We will 
compute the force between the polarized particles by the multiple image method |14]]. We 
will show that the induced force is inversely proportional to the square of the product of the 
relaxation time and the angular velocity of the rotating particle. 
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II. DIPOLE MOMENT OF A ROTATING DIELECTRIC SPHERE 



In this section, we consider two spherical particles, one of which is held fixed, and the 
other one rotates around the axis perpendicular to the line joining the particles centers 
(Fig.0). Both particles are dielectric spheres of dielectric constant ei, suspended in a host 
medium of dielectric constant e m . The rotating (fixed) sphere has a radius a (b). In this 
section, we neglect the mutual polarization between the particles. This assumption is valid 
if the distance between the particles is large. Upon the application of an electric field, 
polarization charges are induced on the surface of the spheres and there are induced dipole 
moments. Suppose the applied field is given by E =E z, the induced dipole moment is given 
by Pao=Paoz, where 

e l e m \ 3 



p a0 = e m E — — J a , (1) 

and a similar expression for the fixed sphere. In the presence of a rotational motion, a 
displacement of polarization charge occurs on the particle surface. 

Suppose the sphere is undergoing an anti-clockwise rotation with an angular velocity 
uo = —ujy. The induced dipole moment will rotate with the same angular velocity, and the 
rate of change of the induced dipole moment is given by: 

§=3xjT (2) 

However, the induced dipole moment can relax back to its original orientation (along Eq) 
and magnitude. If the relaxation process is characterized by a relaxation time r r , we have: 

When a steady state is reached, we have: 

uj x p- —(p- p a0 ) = 0. (4) 
TV- 
Solving for the above equation, we obtain: 
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V 



PaO + T r (u X p a0 ) 



(5) 



1 + (iOT r ) 2 

The result of Eq.fl^) has a simple geometric interpretation. It can be shown that the angle 
between p a0 and p is equal to (ft and the magnitude of p is given by \p\ = p a o cos 0, where 
is given by tan = ujr r . 



III. MULTIPLE IMAGE DIPOLES FOR A PAIR OF DIELECTRIC SPHERES 



In this section, we generalize the multiple image method |L4] to calculate the force 
between the polarized spheres, when the induced dipole moments are not parallel. Consider 
two dielectric spheres, the radii of which are given by a and b respectively and the separation 
between their center is r. Suppose both the spheres are polarized with induced dipole 
moments p a o and p^ respectively. Due to mutual polarization effect, the dipole moment of 
each sphere will induce an image dipole on the other sphere. The induced image dipole will 
further induce another image dipole on the original sphere and hence multiple images are 



formed |fL4|| . To express the image dipole moment, it is convenient to separate the original 
dipole moments into x (perpendicular to E ) and z components (parallel to E Q ): 



PaO = PaOxX + PaOz?, PbO = PbOxX + Pb 0z Z. 



(6) 



The total dipole moments are given by: 



Pa = PaxX + p az Z, Pb = Pbx% + PbzZ- 



(7) 



The x-components can be found by using the multiple image method |L4 

p ate .(&sinha) 3 (-r) 2n_2 PbOx{a sinh a) 3 (-r 
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The ^-components can be found similarly [ITj]: 
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(10) 
(11) 



(asinhna + 6sinh(n — l)a) 3 ' (rsinhna) 3 

The factor (— r) and (2r) apply to the transverse and longitudinal image dipole moments 
respectively. The factor r is known as the dielectric contrast and is given by: 



The parameter a is given by: 



cosh a 



£i ~ e r 
ei + e r 



r 2 — a 2 — b 2 
2ab 



(12) 



(13) 



We should remark that the present generalization is an approximation only, because there is 
a more complicated image method for a dielectric sphere fll5| . However, in the limit r — > 1, 
the above expressions reduce to the results of perfectly conducting spheres. We expect that 
this approximation will be good at high contrast (r — > 1). We have checked the validity of 
the approximation by comparing the above analytic expressions with the numerical solution 



of an integral equation method Hlg] . The force between the spheres is given by the following 



expression, based on the energy considerations [pT7|.p~S 

-t 1 



-V 



E -(p a + p b ) =F r f + F e e. 



(14) 



It should be remarked that F(9) = F{9 + 7r) by examining Fig.|l|. We will calculate the 
results for an anti-clockwise rotation (ou = —uy) only. The results for a clockwise rotation 
can be found similarly. 



IV. NUMERICAL RESULTS 

In this section, we report the separation and angular dependence of the force between 
the spheres. We will calculate the radial component (F r ) and the tangential component (F$) 
of the force and discuss their dependence on ur r . Moreover, we let b = a for convenience. 



In Fig.|2|, we plot the radial force between the spheres against the separation parameter 
a, defined by a = r/(a + b). Two cases are considered, namely, the transverse and longi- 
tudinal field cases. In the transverse (longitudinal) field case, the line joining the spheres 
is perpendicular (parallel) to the applied field. We first examine the transverse field case, 
with uT r being chosen to be 0.1, 1 and 10 respectively. The dielectric contrast is chosen to 
be r = 1/3 (low contrast) and r = 9/11 (high contrast). From Fig.§, our results show that 
the rotational motion of the particle generally reduces the force between the particles. In 
order to see the effects more clearly, we plot in Fig|| the products Ft<J A and F^cr 4 against 
a. It can be seen that for a large separation (a > 3), these quantities tend to constant, 
indicating that the force varies as a -4 . For a small separation (cr < 1.5), the magnitude of 
the transverse force increases rapidly. The longitudinal field case shows a similar behavior: 
the rotational motion generally reduces the magnitude of the interparticle force. 

We next discuss the angular dependence of the interparticle force (Fig.f|). We concentrate 
on the radial component (F r ) first. We can see that, for a small angular velocity (uT r = 0.1), 
the radial component attains its minimum and maximum value at 9 pa and 9 pa tt/2 
respectively, which correspond to the longitudinal and transverse field case (FigfJ). For the 
tangential force, Fq pa at 9 = and 9 = tt/2. This is expected because when the angular 
velocity of the particle is small, the effect of rotation can be neglected. However, when the 
sphere rotates faster {ujt t = 1), the peak shifts towards the direction of smaller 9. Similar 
conclusion can be drawn for the tangential force. In the lower panels of Fig.^, the tangential 
force no longer vanishes at 9 = and 9 = tt/2 for ujt t = 1. From the above results, we 
conclude that the nonrotating particle is forced to move in the direction of increasing 9. 
More precisely, if originally the two particles are at 9 = n/2 while there is no rotational 
motion, the particles will repel each other and there is no tangential motion. However, if 
one of the particle begins to rotate, a tangential force will act on the other particle, forcing 
it to move downwards, since Fq > at 9 = ir/2. 

In the opposite limit, when the angular velocity is large (ur r = 10), the peak of the 
radial force will return to the original position (9 = and 9 = 7r/2), and Fq = when 9 = 

6 



or 9 = tt/2. Furthermore, the magnitude of the force is reduced to nearly one half of its 
original value (Fig.[|). It is because when the angular velocity is large, from p' a0 = p a0 cos(f) 
and tan</> = ujt t (Figfj]), the magnitude of the steady-state dipole moment will be small. 
The situation is reminiscent of a polarized sphere interacting with a dielectric sphere with 
a vanishing polarization charge. 

V. DISCUSSION AND CONCLUSION 

Here a few comments on our results are in order. In this work, only one of the particles 
has been assumed to rotate. In a realistic situation, all the particles can rotate. As our 
theory is general, it can readily be applied to the case when all the particles are in rotational 
motion. 

So far, our proposed relaxation time has no microscopic origin. If the relaxation process 
is originated from a finite conductivity of the particle or host medium, then we can estimate 
the relaxation time based on the Maxwell- Wagner theory of leaky dielectrics [TjJ. For a 
(nonrotating) spherical inclusion embedded in a host medium, the expression is 

r r = (ex + 2e m )/ Oi + 2a m ), 

where a±, a m denote the conductivity of the sphere and host medium respectively. 

We may extend the Maxwell- Wagner theory to polarization relaxation of rotating par- 
ticles. In this case, we should add a term ppv to the polarization current density, where 
pp is the polarization charge density and v = uj x r is the rotating velocity. However, it 
is not possible to convert the extra term into a dielectric constant and the generalization 
becomes more complicated due to the nonuniform polarization charge density inside the ro- 
tating spherical inclusions. We are currently examining the solution of the more complicated 
boundary-value problem. 

However, we can make a crude estimate of the relaxation time. According to the geometry 
of the rotating sphere, the relaxation of polarization charge can be along the particle-fluid 
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interface, so 

l/r r = C r 1 /e 1 + a m /e m . 

When ei is large compared with e m and a m = (insulating host medium), both the expres- 
sions for r r reduce to the same result. 

In conclusion, we have developed a model to calculate the interparticle force between a 
pair of dielectric spheres, in which one sphere is held fixed and the other is in rotational 
motion. We have calculated the force between the spheres by using the multiple image 
method. We have found that the effect of rotation is large when the angular velocity of 
the sphere is large {ojt t » 1). To investigate the dynamic ER effects, we have solved a 
differential equation which describes the dynamic polarization of the rotating particles. The 
theoretical problem described here can be realized in experiments ||13||. 
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FIGURES 

FIG. 1. A rotating polarized dielectric sphere interacting with another non-rotating polarized 
dielectric sphere. 

FIG. 2. Interparticle force Ft and Fl plotted against a. The subscript T denote the transverse 
dipole (9 = 7r/2) and L the longitudinal dipole (6 = 0). The rotational motion of the particle 
generally reduces the interparticle force. 

FIG. 3. Similar to FigJ|, but for F T a 4 and F L a A . The magnitude of the forces increases rapidly 
when the separation becomes small. 

FIG. 4. The angular dependence of force. The angle is expressed in unit of tt, and the separation 
parameter is taken to be <r=l.l. The peaks shift to the direction of smaller 9. For lot t = 1, the 
angular locations 8 = and 9 = n/2 are no longer the equilibrium position of the tangential force. 

FIG. 5. The angular dependence of the magnitude of the force. The magnitude of the force is 
reduced to nearly one half of the original value. 
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